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The Laser Interferometer Space Antenna (LISA) will detect thousands of gravitational wave 
sources. Many of these sources will be overlapping in the sense that their signals will have a non- 
zero cross-correlation. Such overlaps lead to source confusion, which adversely affects how well we 
can extract information about the individual sources. Here we study how source confusion impacts 
parameter estimation for galactic compact binaries, with emphasis on the effects of the number of 
overlaping sources, the time of observation, the gravitational wave frequencies of the sources, and 
the degree of the signal correlations. Our main findings are that the parameter resolution decays 
exponentially with the number of overlapping sources, and super-exponentially with the degree of 
cross-correlation. We also find that an extended mission lifetime is key to disentangling the source 
confusion as the parameter resolution for overlapping sources improves much faster than the usual 
square root of the observation time. 



I. INTRODUCTION 

It is anticipated that the Laser Interferometer Space 
Antenna (LISA) will start collecting data in the next 
decade. LISA will operate in the source-rich low fre- 
quency portion of the gravitational wave spectrum, where 
the primary sources are expected to be compact galac- 
tic binaries, supermassive black hole binaries, and ex- 
treme mass ratio inspirals (such as a white dwarf falling 
into a supermassive black hole). A natural question to 
ask is how well the source parameters, such as sky lo- 
cation, masses and distances, can be recovered from the 
LISA data stream. Estimates of the parameter resolution 
have been given for compact galactic binaries 0, 0> @, 
supermassive black hole binaries 0, El and ex- 
treme mass ratio inspirals [9]. These studies focused on 
the problem of identifying one source at a time, and did 
not address the problem of source confusion. The large 
signal from the galactic population of close white dwarf 
binaries was not ignored, but it was treated as an addi- 
tional source of stationary, Gaussian noise to be added 
to the instrument noise. We will see that this is not a 
good approximation. 

Here we study how parameter estimation is affected by 
the presence of multiple overlapping signals in the LISA 
data stream. Importantly, we find that the behavior is 
very different from what one would predict by treating 
the additional sources as stationary, Gaussian noise. The 
problem of source confusion will be most pronounced be- 
low ~ 2 niHz, where it is expected that many tens to 
tens of thousands of galactic binaries will have signals 
that overlap in each frequency bin 0]. Thus, we will fo- 
cus our attention on low frequency galactic binaries that 
are close to one another in frequency. 

The paper is organized as follows. We begin with a 
summary of our findings in Section^ This is followed by 
a review of parameter estimation for an isolated galactic 
binary in Section IlIII Parameter estimation with multi- 
ple sources is described in Section llVI and our plan for 
exploring the parameter space is described in Section Ivl 
The main results are described in Section IVTI A brief in- 



formation theory perspective is presented in Section [VIII 
and concluding remarks are made in Section IVIIII 

II. SUMMARY OF RESULTS 

We find that the parameter estimation uncertainties 
grow exponentially with the number, N, of overlapping 
sources. This should be contrasted to the ViV increase 
one would predict if the other sources were treated as sta- 
tionary, Gaussian noise. The degredation in resolution 
was found to be nearly uniform across the seven parame- 
ters that LISA will measure for a galactic binary system. 
As one might expect, the parameter uncertainties are a 
strong function of the signal cross-correlations. We find 
that as the two signals become more correlated, the pa- 
rameter estimation uncertainties grow at a rate that is 
faster than exponential. Importantly, we find that the 
parameter uncertainties decrease rapidly as the observa- 
tion time T is increased - far more rapidly than the usual 
1/Vr improvement one expects when competing with 
stationary, Gaussian noise. 

III. REVIEW OF PARAMETER ESTIMATION 

Typical galactic binaries can be treated as circular and 
monochromatic. They are described by seven parame- 
ters: sky location (#,</>); gravitational wave frequency /; 
amplitude A; inclination and polarization angles {i,ip); 
and the initial orbital phase 7. The response of the LISA 
instrument to a gravitational wave source is encoded in 
the (Micheleson-like) X(t) and Y(t) time-delay interfer- 
ometry variables [13 ■ We employ the rigid adiabatic ap- 
proximation 01 to describe these variables, and work 
with the orthogonal combinations 3] 

S I (t)=X(t), S n (t) = -L(X(t) + 2Y(t)) . (1) 

Each signal is a function of the parameters A — > 
(8, <fi, /, A, 1, ip, 7) that describe the source. Denoting the 



2 



noise spectral power by S n (f), we adopt the usual noise- 
weighted inner product 



i\b) = 2 



a*(f)b(f)+a(f)b*(f) 

Sn(f) 



df. (2) 



Decomposing the output channels into signal h and 
noise n, the signal to noise ratio is given by 



SNR 2 = (MM- 



(3) 



a=I,II 



For large SNR the parameter estimation uncertainties 
AA ! will have the Gaussian probability distribution ^2] 



det(r) 



2tt 



■ exp 



--I^-AA'AA' ) , (4) 



where the Fisher information matrix T is defined by 



(5) 



=i,ii 



Here h t i = dh/d\ l . For large SNR the variance- 
covariance matrix is given by 



(0) 



and the uncertainties in the parameters are given by 
A A 1 = (C") 1 / 2 . The volume of the n — a uncertainty 
ellipsoid in the d dimensional parameter space is 



TT d / 2 n d , 



r(d/2 + 1) 



(7) 



Applying the above formalism to an isolated galactic 
binary yields the parameter uncertainties shown in Fig- 
ures Q] and |3 The plots were generated by taking the 
median values for 10 5 different sources, each normalized 
to a signal-to-noise ratio of SNR = 10. The uncertainty 
in the frequency has been increased by a factor of 10 6 for 
plotting purposes. 

As can be seen in Figure ^ the uncertainties in five 
of the seven parameters are fairly constant across the 
LISA band. The two exceptions are the sky location 
variables (0,(f>), which show a marked decrease in their 
uncertainties above 1 mHz. This behavior can be traced 
to the time varying Doppler shift, which is one of the 
two ways LISA is able to locate sources on the sky. The 
Doppler shift increases linearly with /, which translates 
into a 1/ f decrease in the positional uncertainties above 
1 mHz. Below this frequency the the angular resolu- 
tion comes mainly from the time varying antenna sweep, 
which is weakly dependent on / below the transfer fre- 
quency /» ~ 10 mHz. 
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FIG. 1: Median parameter uncertainties versus frequency for 
isolated monochromatic binary sources with one year of ob- 
servation and fixed signal-to-noise ratios of SNR = 10. The 
parameters are /, In A, 9, <j>, t, ip, 7. 
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FIG. 2: Parameter uncertainties versus time for isolated 
monochromatic binary sources with f — 1 mHz and fixed 
signal-to-noise ratios of SNR = 10. The parameters are 
In /, In A, 6, <f>,i, V>,7- 



As can be seen in Figure [21 the uncertainties in six 
of the seven parameters decrease as 1/ vTofcl for obser- 
vation times longer than a year, in accordance with the 
standard \JT h s increase in the SNR. The one exception 
is the frequency, which benefits from an additional 1/T b s 
shrinkage in the width of the frequency bins, leading to an 
overall decay in the frequency uncertainty. The 

shrinkage in the size of the frequency bins is illustrated in 
Figure |3 where the discrete Fourier transform of a typ- 
ical signal is shown after one year and after ten years of 
observation. Later we will see how this improved resolu- 
tion of the sidebands helps reduce source confusion when 
multiple overlapping sources are present. For observa- 
tion times less than one year the positional uncertainties 

also decay as l/T^ 2 . The additional factor of l/T a b s 
can be traced to the two mechanisms by which LISA can 
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FIG. 3: Fourier transforms of a typical signal for observation 
times of one and ten years. 



locate a source: antenna sweep and Doppler shift. The 
angular resolution due to antenna sweep depends on the 
angle through which the antenna has been swept, which 
increases as T i, s for T obs < year. The angular resolution 
due to the variable Doppler shift can be thought of as a 
form of long baseline interferometery. The resolution im- 
proves with the length of the synthesized baseline, which 
grows as T obs for T obs < year. 

In what follows, we will typically quote our results 
in terms of uncertainty ratios. These ratios will com- 
pare the parameter resolution when one or more overlap- 
ping sources are present to the parameter resolution that 
would be possible if each source were isolated. To arrive 
at the absolute uncertainties one needs to multiply by 
the isolated source uncertainties quoted above. 



IV. MULTIPLE SOURCES 

The Fisher information matrix approach to estimating 
parameter uncertainties is easily generalized to multiple 
sources. For TV circular, monochromatic binary systems, 



the parameter space is 7N dimensional, and we have the 
parameter vector 



A = A 1 + A 2 + . . . + A 



N 



(8) 



We will use the notation A A to refer to all 7 parameters 
of the A th binary. The total signal S — Si or S = Sn 
will be the sum of the signals for the individual binary 
systems, 



S = S 1 + S 2 + S 3 



S 



N 



(9) 



Here S A is the signal due to the A th binary. Since 
dS A /dA B = for A ^ B, the full Fisher information 
matrix has a simple structure. If we define 
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the diagonal blocks V aa ar e the usual 7x7 Fisher infor- 
mation matrix for the A th source, while the off-diagonal 
blocks Tab describe how the parameter estimation for 
source A is influenced by the presence of source B and 
vice- versa. If all the sources are un-correlated (for exam- 
ple, they might be widely spaced in frequency), then the 
off-diagonal blocks will be zero, and the full Fisher in- 
formation matrix will be block diagonal. Upon inverting 
to get the variance-covariance matrix the block diagonal 
structure will be preserved, with each 7x7 block equal 
to the inverse of the corresponding single source Fisher 
information matrix. However, if the sources are over- 
lapping, the off-diagonal blocks will be non-zero, which 
will affect the values of the diagonal elements in the full 
variance-covariance matrix. The number of off-diagonal 
elements grows quadratically with the number of sources, 
while the number of diagonal elements grows linearly 
with the number of sources. Thus, we anticipate that 
the ability to resolve a particular source's parameters will 
degrade rapidly as the number of correlated sources in- 
creases. 

We also anticipate that the parameter uncertainties 
will depend strongly on the degree of correlation between 
sources. The correlation between two signals S A and S B 
is defined: 



(S A \S B ^ 
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(11) 



In analyzing the results, two similar approaches will be 
used to obtain a quantitative measure of the increase in 
the parameter estimation uncertainties. First is a global 
comparison of the uncertainties. The uncertainties form 
an ellipsoid in the parameter space whose volume is given 
by the determinant of the variance-covariance matrix. 
One measure of the uncertainty increase due to the cor- 
relation of the binary systems is given by a ratio of the 
geometric mean of the uncertainties (GMUR): 
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GMUR = 



dctr 



1/7N 



(12) 



The GMUR describes the mean increase in the parameter 
uncertainties due to source confusion. 

A second measure of the uncertainty increase is a pa- 
rameter by parameter comparison, looking at the ratio 
of the parameter uncertainty (PUR) between the uncer- 
tainty of a particular parameter calculated in the 7N- 
dimensional variance-covariance matrix, AA^, and the 
uncertainty of the same parameter calculated in the iso- 
lated binary's variance-covariance matrix, AA^. 
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(13) 



This definition for the PUR is independent of the am- 
plitudes of the sources, as can be seen in the simple case 
of two sources, each described by one a. The PUR for 
ai is given by: 



PUR 



AA r 



AA ai 



V 2 
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(14) 



where 



^/Paiai Ta2Q2 



(15) 



The independence of the degree of confusion on the 
source amplitudes would seem to imply that an arbitar- 
ily weak source can affect parameter estimation to the 
same degree as an arbitarily strong source. Clearly, this 
makes no sense in the limit that the amplitude of the 
weak source goes to zero. The resolution to this apparent 
paradox is that the Fisher information matrix approach 
is only meaningful for sources with SNR > 1, so arbitar- 
ily weak sources are not permitted in the analysis. 



V. SOURCE PARAMETER SELECTION 

The parameter space for A slowly evolving binaries is 
7 A dimensional. This large dimensionality makes it dif- 
ficult to carry out an exhaustive exploration of all the 
circumstances that can affect parameter estimation. We 
focused on sources that were close in frequency, and chose 
the other source parameters randomly. The frequency 
chosen for the first binary fixes the base frequency /base 
of a data run. The frequencies of the remaining A — 1 bi- 
naries were assigned frequencies of /basc+(*+^)/m, where 
fm = 1/year is the modulation frequency, i is an integer, 
and a; is a random number between and 1. Thus the 
remaining binaries are between i and i + 1 modulation 
frequency units from the base frequency. Sky locations 
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FIG. 4: Median PUR versus time of observation for two, 
three, and four binaries with base frequencies of 1 mHz from 
an all sky draw. 



were chosen by two methods, the first being a random 
draw on cos(#) G [—1, 1] and <fi S [0, 2ir) with a fixed dis- 
tance to the binary of 1 kiloparsec, the second being a 
random draw from a galactic distribution of 9, (j), and 
binary distance. The values of the polarization ip and or- 
bital phase 7 were randomly drawn between and tt and 
and 2tt respectively. The inclination l was taken from a 
random draw on cos(t) G [—1,1], but with values outside 
of the range l £ [1°, 179°] rejected to avoid the degen- 
eracy in the gravitational wave produced by a circular 
binary viewed along its axis. Masses were taken to equal 
0.5Mq, as we wished to model white dwarf binaries. 



VI. RESULTS 
A. Varying time of observation 

The current lifespan of LISA is estimated to be five 
years. Our data show that the longer LISA observes the 
binaries systems, the less the effect of confusion between 
sources. Figure 0] shows the median parameter uncer- 
tainty ratios (PURs) for two, three, and four binary sys- 
tems with base frequencies of 1 mHz as a function of 
the time of observation. Figure shows the median ge- 
ometric mean uncertainty ratio (GMUR) for the same 
data. Each set of PUR data in Figure 0] contains the un- 
certainty ratios for all seven parameters, and as can be 
seen in the figures, the uncertainty ratios for each of the 
seven parameters are closely related. This relationship 
holds true for other frequencies in the LISA band. It is 
this uniformity in the increase in parameter uncertainties 
that motivates our use of the GMUR as an estimate for 
the individual PURs. Thus, subsequent results will be 
given in terms of the GMUR. 

Figure shows how the median GMUR increases with 
the number of binaries for one year of observation, at 
base frequencies of 0.1, 1, and 5 mHz. Figure shows 
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FIG. 5: Median GMUR versus time of observation for two, 
three, and four binaries with base frequencies of 1 mHz from 
an all sky draw. 
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FIG. 7: Median GMUR plotted against the number of binary 
systems for one, two, and three years of observation at a base 
frequency of 1 mHz. 
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Number of Binary Systems FIG. 8: Median Correlation magnitude versus time of obser- 

vation for two binaries with a base frequency of 1 mHz from 
FIG. 6: Median GMUR plotted against the number of binary a n all sky draw, 
systems for one year of observation at base frequencies of 0.1, 
1, and 5 mHz. 

B. Uncertainty ratios for two binary systems as a 
function of frequency difference 



how extending the time of observation affects the GMUR. 
In both instances, the parameter uncertainties increase 
roughly exponentially with the number of overlapping 
sources. 

Figurc|Hlshows the median of the magnitude of the cor- 
relation between the signals of two binary systems with 
a base frequency of 1 mHz taken from an all sky draw. 
As expected, the correlation magnitude falls off as l/T Q b s 
for observation times greater than a year (the numerator 
in ltTT|) oscillates, while the denominator grows as T b s )- 
The 1 /T b s fall-off in the signal correlation should be con- 
trasted with the much faster fall-off in the GMUR seen 
in Figure The reason for this difference in fall-off will 
become clear in Section IVl CI 



Figure shows how the median parameter uncertain- 
ties depend on the frequency separation when two binary 
systems are present. The base frequency, /& ase , for the 
first binary was held fixed at 1 mHz, while the frequency 
of the second binary was randomly chosen between fbase 
and fbase+ifm for i between and 20. Figure|Hlshows the 
plots for GMUR versus i for all sky and galactic draws 
for one year of observation. 

The uncertainty ratio drops rapidly as the frequency 
separation increases, approaching unity by the time the 
binaries are ~ 5/ m apart. The frequency difference of ~ 
5/ m corresponds to the typical half-bandwidth of a source 
at 1 mHz. In other words, 1 mHz sources seperated by 
> 5/ m have almost no overlap, and thus there is little 
source confusion for i > 5. 
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FIG. 9: Median GMUR plotted against modulation frequency 
bin difference of two binaries from an all sky and a galactic 
draw with a base frequency of 1 mHz and one year of obser- 
vation 
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FIG. 10: The median GMUR plotted against signal corre- 
lation of two binaries from a galactic sky draw with a base 
frequency of 1 mHz and one year of observation 



C. Uncertainty ratios as a function of signal 
correlation 

Figure ITU1 plots the median GMUR versus the signal 
correlation for a galactic draw when two binaries are 
present within Sf = f m of each other. The plot shows 
that the GMURs increase faster than exponentially as the 
signal correlation increases. While a link between signal 
correlation and uncertainty ratios is to be expected, the 
super-exponential nature of the relationship came as a 
suprise. Figures ITU and IT21 show the spreads for correla- 
tion and GMUR in the same data run. As can be seen 
in the log-linear scale of Figure ^1 the GMUR falls off 
exponentially from its most frequent value of 1.2. 

While Figure ITT1 shows the results for two correlating 
binaries, its results can be extended to multiple correlat- 
ing binaries. For example, in the figure one can see that 
87% of the pairs of binaries have correlation magnitudes 
|k| < 0.5. Thus with N binaries the probability that 
there will be a pair of binaries with |k| > 0.5 grows as 
1 — o.gjiNiN- 1 )/ 2 ). As binaries overlapping in frequency 
are expected to number from tens to tens of thousands, 
the probability of having high correlation is very close to 
1. 

The rapid increase in the parameter uncertainties 
with signal correlation explains the results seen in Sec- 
tion IVI Al The 1/T f, s decrease in the signal correlation 
translates into a much faster decrease in the GMUR as a 
function of observation time. 



D. Uncertainty ratios as a function of base 
frequency 

Figure ITTjl shows how the parameter uncertainty ratios 
depend on base frequency. The plots are for all sky and 
galactic draws with two binaries for one year of obser- 
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FIG. 11: A histogram the correlation of two binaries from a 
galactic draw with a base frequency of 1 mHz and one year 
of observation 
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two binaries from all sky and galactic draws with one year of 
observation 
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FIG. 14: Median GMUR plotted against azimuthal sky sepa- 
ration for two pairs of binaries from a galactic draw with one 
year of observation with base frequencies of 0.1 and 1 mHz. 



vation. The base frequency, fbase, was varied between 
0.01 and 10 mHz, while the frequency of the second bi- 
nary was randomly chosen between fbase and fbase + 
f m . The GMUR was found to be fairly constant be- 
low 1 mHz, followed by a rapid decrease at frequencies 
above 1 mHz. This behavior can be traced to the different 
Doppler shifts experienced by each source. The motion 
of LISA relative to a source located at (9, <fi) imparts a 
Doppler shift equal to 



Sf D ~ 7rsin0sin(27r/ m t - <j>) [j^^jf™- ( 16 ) 

The magnitude of 8fn becomes comparable to the fre- 
quency resolution Sf = l/T b s = 1/year for / ~ 0.3 
mHz. Sources that are well seperated in ecliptic azimuth 
<j) will experience Doppler shifts that differ in sign, as 
LISA will be moving toward one source and away from 
the other. Thus, we expect that the degree of source con- 
fusion will depend on the azimuthal separation of the two 
sources. Our expectations are confirmed in Figure 1141 
where the dependence of the GMUR on azimuthal sepa- 
ration is plotted for base frequencies of 0.1 and 1 mHz. 
Unfortunately, most galactic sources are within ~ 20° 
degrees of each other in ecliptic azimuth, which helps ex- 
plain the larger GMURs for the galactic distribution as 
compared to the all-sky distribution. 



VII. INFORMATION THEORY PERSPECTIVE 

There have been several attempts to understand LISA 
source confusion in the framework of information the- 
ory [LH [13 According to Shannon 16], the maxi- 
mum amount of information that can be transmitted over 
a noisy channel with bandwidth B in a time T b s is 



I s =BT ohs \og 2 (l+^\. (17) 

Here Ph and P n are, respectively, the signal and noise 
power across the bandwidth. (Note: The SNR ratio is re- 
lated to these quantites by SNR 2 ~ B T^P^/ P n ). As an 
example, a typical bright galactic source at 1 mHz with 
Ph/Pn ~ 10 can transmit at most 7s — 70 bits of infor- 
mation in one year via the two independent LISA data 
channels Si and Sn. In practice the actual information 
content will be somewhat less due to sub-optimal encod- 
ing. It takes log 2 (a;/Aa;) bits of information to store a 
number x to accuracy Ax, thus the total amount of in- 
formation required to describe a monochromatic binary 
to a precision AA is 




Using the results in Figure^for a bright galactic source 
at 1 mHz with Ph/P n ~ 10, we find that 7a — 45, which 
indicates that the encoding, while sub-optimal, is quiet 
respectable. 

One might hope that information theory could be used 
to predict some of the other results described in Sec- 
tion Unfortunately, this proves not to be the case, 
for while information theory sets limits on how well one 
might do, it sets no limits on how poorly. The scaling 
of the GMUR as a function of the number of sources is 
a good example. The information needed to localize N 
galactic sources is equal to 
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V 7 



N 



7JV 



(19) 



VJ* is the volume 



where AVtn is given by J7| and Vjn — y? 
of the 7N dimensional parameter space. Requiring that 
Ia{N) < Is(N) yields a lower bound on the geometric 
mean of the parameter uncertainties (GMU) of 



GMU > 



(1 + P s /P n 



BT/7N 



(20) 



Consider first the case of N isolated binaries of similar 
brightness. The power ratio P s /P n will be independent 
of N, and the bandwidth will be the sum of the band- 
widths of each signal, so B will grow linearly with N. 
Thus, equation (121 H implies that for isolated binaries, 
the parameter uncertainties will be independent of the 
number of sources. This is indeed the case, as the Fisher 
information matrix is block diagonal for isolated sources. 
Now consider the case of N overlapping binaries sharing 
the same fixed bandwidth B. The power P s will grow 
approximately linearly with N , leading to the prediction 
that the GMU must grow faster than (l/iV) (1/Ar) . The 
actual scaling seen in Figures[£]and[7|tcll a different story, 
with the GMUR increasing as ~ N N . Thus, while infor- 
mation theory provides a lower bound on how quickly the 
parameter uncertainties must increase with the number 
of overlapping sources, the bound is too weak to be of any 
real use. The weakness of the bound is probably related 



to how poorly the information for multiple overlapping 
binaries is encoded. 



VIII. DISCUSSION 

We have found that source confusion will significantly 
impact LISA's ability to resolve individual galactic bina- 
ries. We found that source confusion affects parameter 
estimation fairly uniformly across the seven paramters 
that describe a galactic binary. 

The two most significant findings were that source con- 
fusion grows exponentially with the number of correlated 
sources, and that source confusion decreases rapidly with 
time of observation. 

Source confusion will be a significant problem in the 
frequency range between 0.01 and 3 mHz, where it is 
estimated that there are upwards of 10 8 galactic binary 
systems. The decrease in the parameter uncertainties 
with time of observation is far greater than the usual 
1/y/Tobs decay that occurs when competing against sta- 
tionary, Gaussian noise. The fast improvement is due to 
the 1 /T }, s decrease in the source cross-correlation, and 
the sensitive dependence of the parameter uncertainties 
on the degree of cross-correlation. 

Our findings suggest that earlier work that treated the 
galactic background as stationary, Gaussian noise should 
be revisited, and that every effort should be made to 
extend the LISA mission lifetime beyond its nominal 3 
year duration. 
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